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Abstract 

We use chiral effective theory (xET) to predict the deuteron form factor ratio 
Gc/Gq as well as ratios of deuteron to nucleon form factors. These ratios are 
calculated to next-to-next-to-leading order. At this order the chiral expansion for 
the NN isoscalar charge operator (including consistently calculated 1/M correc- 
tions) is a parameter-free prediction of the effective theory. Use of this operator in 
conjunction with NLO and NNLO xET wave functions produces results that are 
consistent with extant experimental data for < 0.35 GeV^. These x^T wave 
functions predict a deuteron quadrupole moment Gq{Q^ = 0) = 0.278-0.282 fm^ — 
with the variation arising from short-distance contributions to this quantity. The 
variation is of the same size as the discrepancy between the theoretical result and 
the experimental value. This motivates the renormalization of Gq via a two-nucleon 
operator that couples to quadrupole photons. After that renormalization we obtain 
a robust prediction for the shape of Gc/Gq at < 0.3 GeV^. This allows us to 
make precise, model-independent predictions for the values of this ratio that will 
be measured at the lower end of the kinematic range explored at BLAST. We also 
present results for the ratio Gq/Gm- 
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1 Introduction 



The static electromagnetic properties of deuterium provide interesting information on the 
dynamics at work within the nucleus. The fact that deuterium's charge is one teaches us 
little other than the validity of charge conservation in the nuclear system, but that its 
magnetic moment /id 7^ /in + /ip and that it has of a non-zero quadrupole moment are 
facts which played an important role in establishing that non-central components of the 
A^A^ potential are at work within the deuterium nucleus. 

The most accurate value for the deuteron quadrupole moment comes from a molecular 
physics experiment 0121. It is: 



Qd = 0.2859(3) fml 



(1) 



Meanwhile the best determination of fid comes from the spectroscopy of the deuterium 
atom. It is [3]: 

/irf = 0.8574382284(94) fi^. (2) 

But elastic electron scattering from deuterium — provided the one-photon-exchange ap- 
proximation is valid — probes Ml and E2 responses for (virtual) photons that have a finite 
three-momentum q. The relevant form factors are related to Breit-frame matrix elements 
of the two-nucleon four-current via 
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provide a complete set of invariant functions for the description of the deuterium four- 
current that interacts with the electron's current in this approximation. In Eqs. 
we have labeled the deuteron states by the projection of the deuteron spin along the 
direction of the three- vector Pg — pe = q, and 77 = Q^/(4Mj), with = |qp since we 
are in the Breit frame. We can then calculate the deuteron structure functions: 
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In terms of A and B the one-photon-exchange interaction yields a lab. frame differential 
cross section for electron-deuteron scattering 
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Here 6^ is the electron scattering angle, and ^j^g is the (one-photon-exchange) cross 
section for electron scattering from a point particle of charge |e| and mass M^. The form 
factors defined in Eqs. ©^(0) are related to the static moments of the nucleus by: 

Gc(0) = 1, (9) 
Gq{Q) = Qd, (10) 

Gm{0) = /i.^, (11) 

with M the nucleon mass. For recent reviews of experimental and theoretical work on 
elastic electron-deuteron scattering see Refs. [H El E] • 

From Eq. (jH)) it is already clear that measurements of the differential cross section alone 
cannot yield uncorrelated information on all three form factors. To measure Gq, Gm, and 
Gc in a model-independent way one must obtain data with polarized deuterium targets 
or polarized electron beams. A new set of measurements of polarization observables in 
electron-deuteron scattering will soon be available from the data set obtained at the Bates 
Large- Acceptance Spectrometer Toroid (BLAST). There polarized electrons of energy 850 
MeV circulated in the Bates ring and were scattered from an internal target containing 
polarized deuterium. The significant amount of beam on target (3 million Coulombs since 
late 2003), and high degree of beam and target polarization achieved at BLAST, means 
that we anticipate data on electron-deuteron polarization observables that is more precise 
than that derived from any previous measurement. 

The electron beam circulating in the Bates ring was roughly 70% polarized, and the 
deuterium target employed could operate in both a vector-polarized and tensor-polarized 
mode. This gives access to all of the elastic electron scattering deuterium structure 
functions. Prominent among these are tu, the vector analyzing power, and ^20, the tensor 
analyzing power. They are related to the form factors defined above by 

t - ,/o ^(^ + 2) + |y^(l + 25) . _ 1 + f 

- '"^^TTW+WTW "" ^'^1 + 2^+7(1 + 2^))' ^ ^ 



where the ratios x and y are: 
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and e = {1 + 77) tan 2 j' Therefore measurements of tn and t2o should facilitate the 
extraction of the ratios Gm/Gc (which at the Q^'s we will consider here mainly affects 
til) and Gq/Gc (which mainly affects ^20)- In this paper we provide predictions for these 
ratios which are based on chiral effective theory (^ET). 
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This approach (for reviews see Refs. [71 |H1 IHl HO]) is based on the use of a chiral 
expansion for the physics of the two-nucleon system. In the formulation suggested by 
Weinberg [TT| IT^ IT^ . the xET treatment of the A^A^ system is based on a systematic 
chiral and momentum expansion for the two-nucleon-irreducible kernels of the processes 
of interest. In particular, wave functions are computed using an A^A^ potential expanded 
up to a given order in the small parameter: 



where p is the momentum of the nucleons and A is the breakdown scale of the theory. 
For electron-deuteron scattering the other two-nucleon-irreducible kernel that must be 
calculated is the deuteron current operator J^. We also expand this object as: 



where the operator (9^*^ contains i — 1 powers of the small parameter P, which now includes 
the momentum transfer to the nucleus, q, as one of the small scales in the numerator. 
For chiral effective theories without an explicit Delta degree of freedom A will in general 
be Ma — M, but in electron-deuteron elastic scattering the AA^ intermediate state is not 
allowed and so A will be larger, A > 2 (Ma — M). 

The A^A^ potential has now been computed up to 0{P^) 0{P^) 
and 0{P'^) ^71 Ej. In this paper we will employ wave functions computed using the 
next-to- leading order [NLO=0(P^)], next-to- next-to-leading order [NNLO=0(P^)] and 
N^LO potentials [0(P^)] developed in Ref. jTH|. These potentials are regularized in two 
different ways: first, spectral-function regularization (SFR) at a scale A [1^], is applied 
to the two-pion contributions. Then, after the SFR potential Vi1i'{p,p') is obtained in 
a particular A^A^ partial wave, it is multiplied by a regulator function /, so that the 
Lippmann-Schwinger equation can be straightforwardly solved: 



There have been questions raised as to the consistency of the wave functions computed 
in this way [201 1211 1221 1211 1211 1211 • Partly because of these questions we will, for compar- 
ison, also present results for electron-deuteron matrix elements using the form pfi|l for 
the current operator, and wave functions derived from the Nijm93 [22] or CD-Bonn [2Z| 
potentials, as well as potentials with one-pion exchange at long range and a square well 
and surface delta function of radius R [28,. We stress that such calculations are not chi- 
rally consistent. However, common features of deuteron observables that can be identified 
within calculations that use these different types of wave functions — chiral effective the- 
ory, potential models, and one-pion-exchange tails — should be independent of the details 
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of physics at ranges -C 1/ in deuterium, and so should not be sensitive to any subtleties 
pertaining to the renormalization of the xET. 

The operators O^^^ and the coefficients in Eq. (fT^ are constructed according to 
the counting rules and Lagrangian of heavy-baryon chiral perturbation theory (HBxPT), 
which is reviewed in Ref. Here the results we will present for Gc and Gq include 

all contributions to up to chiral order eP^. This is the next-to- next-to-leading order 
(NNLO) for these quantities. Calculations of electron-deuteron scattering with the NNLO 
xET operator were already considered in Ref. [3T], which improved upon results with the 
O(eP^) operator in Ref. [HSj and the 0(e) results of Ref. p8|. However, as was already 
observed in Ref. [HI] and is reiterated below, calculation of the quadrupole combination 
of matrix elements at NNLO does not reproduce the experimental value of Qd to the 
accuracy one would expect at that order. We identify the cause of this as short- distance 
two-body contributions to Jq of natural size (i.e. with ~ 1) through which quadrupole 
photons induce anL = 0— *L = transition in the S = 1 deuteron state |2H1 121] • We use 
the operator induced by these short-distance contributions to renormalize the deuteron 
quadrupole moment, and hence the form factor Gq. We also provide results for Gm up to 
NLO. Not surprisingly, Gm at NLO proves more sensitive to short-distance physics than 
does the renormalized Gq. 

Throughout this work we will use the factorization of nucleon structure employed in 
Ref. in order to include the effects of finite nucleon size in the calculation. There it 
was shown that the chiral expansion for the ratios: 

-^2- and -^S- and -^rr, (18) 

^(■5) M^) -^(s)' ^ ^ 

with Gg'* and G'^^} the isoscalar single-nucleon electric and magnetic form factors, is better 
behaved than the chiral expansion for Gc-, Gq., and Gm themselves. The ratios (|18|) allow 
us to focus on the ability of the chiral expansion to describe deuteron structure, and we 
will employ the xET results for the ratios in our efforts to predict BLAST'S results for 
polarized electron scattering from a deuterium target. We note that, up to the order 
we work to here, our predictions for Gc/Gq are independent of the manner in which we 
include nucleon structure in the calculation. Our invoking the factorization of nucleon 
structure in the electron-deuteron matrix elements plays no role in our predictions for 
Gc/Gq. 

The chiral perturbation theory for this calculation was laid out in Refs. jSHES, and so 
here we merely summarize the pertinent features of the chiral expansion for the deuteron 
currents in Section|21 However, in doing so we find that we must address the issue of how to 
calculate the corrections to this ratio that have coefficients which are fixed by low-energy 
Lorentz invariance. We deal with this problem in Section IHl by recalling results of Friar, 
Adam et al, and Arenhovel et al, which show that such corrections can be calculated 
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unambiguously, as long as they are included consistently in both the NN potential and 
the current operator. Then, in Section IH we discuss effects in the Jo operator beyond 
O(eP^), and explain how we will estimate their impact on Gc and Gq. In particular, 
we write down an operator that represents the effects of physics at mass scales above 1 
GeV on Gq, and can repair the discrepancy between the experimental value of Qd and 
our predictions for Gq{0). We also discuss how to estimate the remaining uncertainty in 
our results. Then, in SectionElwe present results for Gc, Gq, and the ratio Gc/Gq. We 
show that the shape of Gq can be predicted in a model-independent way for < 0.3 
GeV^, but the uncertainty in the ratio Gc/Gq is sizeable at ~ 0.4 GeV^. Finally, in 
Section El we present results for Gc/Gm, and in Section [7| we summarize and provide an 
outlook. 

2 The deuteron current 

We now discuss the charge and current operators in turn. Such a decomposition is, of 
course, not Lorentz invariant, so here we make this specification in the Breit frame. 

2.1 Deuteron charge 

The vertex from C^^)^ which represents an Aq photon coupling to a point nucleon gives 
the leading-order (LO) contribution to Jo as depicted in Fig. IH^a). At O(eP^) this is 
corrected by insertions in £^jv that generate the nucleon's isoscalar charge radius. This 
gives a result for Jo through O(eP^): 

4'^ = \e\{l-l{rh)Q')+jr''\^) (19) 

with jo^^^^ ■* (q) the "relativistic" corrections to the single- nucleon charge operator. These 
contributions have fixed coefficients that are determined by the requirements of Poincare 
invariance. Since these coefficients scale as this particular set of O(eP^) contri- 

butions are generally smaller than one would estimate given the formula ()15p for the 
parameter P. These "relativistic" corrections can be calculated by writing down a Jo op- 
erator that, when inserted between deuteron wave functions calculated in the two- nucleon 
center-of-mass frame, yields results for the matrix elements that are Lorentz covariant 
up to the order to which we work. To do this we employ the formalism of Adam and 
Arenhovel, as described in Ref. |33j . 

Meanwhile the only contribution at 0{eP^), or NNLO, comes from the tree-level two- 
body graph shown in Fig. ^b). In HBxPT the relevant single- nucleon photo-pion vertex 
arises as a consequence of the Foldy-Wouthuysen transformation which generates a term 
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Figure 1: Diagrams representing the leading contribution to the deuteron charge operator 
[(a)], the leading two-body contribution to Jq [(b)], and the dominant short- distance piece 
[(c)]. Solid circles are vertices from C^^j^, and the shaded circle is the vertex from C^^^^- 



in C'^'^\ Straightforward application of the Feynman rules for the relevant pieces of the 
HBxPT Lagrangian gives the result for this piece of the deuteron current p?T] : 



(pVi'Hq)|p) = « 



ai ■ qq-2 ■ (p - p' + q/2) 
m2 + (p - p' + q/2)2 



+ (1^2) 



(20) 



where p and p' are the (Breit-frame) relative momenta of the two nucleons in the initial 
and final-state respectively (see Ref. [33] for a much earlier derivation). 

Thus we now have a result for the deuteron's charge operator which can be summarized 



as: 



(pVo(q)|p) 



\e\[l-l{rl.)Q')+jr'^i 



Si^){p'-p-q/2) + {p'\4'\q)\p)+0{eP 



(21) 

However, it was shown in Ref. that the parameterization (|19p of the nucleon's isoscalar 
charge distribution breaks down at |q| ~ 300 MeV. In order to avoid this difficulty we 
observe that the result ()21|) may be recast, up to the order to which we work, as: 



(pVo(q)|p) 



P 



q/2) + (pV^'^(q)lp)] G^^\q') + 0{eP'), 

(22) 

with G^^\Q'^) the complete one-loop HBxPT result for the nucleon's isoscalar electric form 
factor This means that we can write a result that is independent of xPT's difficulties 
in describing nucleon structure if we focus on the ratio of (p'|^o(q)|p) to G^^\Q'^). We 
will then use experimental data ^, in particular the parameterization of Mergell et al. [37J, 
for g'^^{Q'^) in all computations we present below. This use of experimental data for the 
single-nucleon matrix element that appears in Eq. allows us to focus on how well 
the xET is describing deuteron structure, since it removes the nucleon-structure issues 

^There is, of course, some circularity here, since electron-deuteron elastic scattering data has been 
used to constrain the neutron electric form factor, see, in particular, Ref. |36| . 
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from the computation of the deuterium matrix element. Our technique to achieve this 
is rigorous, up to the chiral order to which we work here. The ratio Gc/Gq can also be 
computed independent of nucleon-structure issues, as is made clear by a brief examination 
of Eq. , together with the definitions (0) and (jH) . 

2.2 Deuteron three-current 

The counting for the isoscalar three-vector current J was already considered in detail by 
Park and collaborators [38j. J begins at 0(eP), but at 0{eP^) there are finite-size and 
relativistic corrections, which are suppressed by two powers of P^. This is the highest 
order we will calculate Gm to here, and at this order, using factorization we have: 

(p'l J(=^)|p) = [|e|(p + q/2)/M + ^^^s<y x q + j(VM^)(q)]Gl? (Q2)5(3)(p' - p - g/2). (23) 

with p — q/4 is the momentum of the struck nucleon, and /i^ is the isoscalar magnetic 
moment of the nucleon, whose value we take to be = 0.88|e|/(2M). 

At O(eP^) [NNLO] two kinds of magnetic two-body current enter the calculation. 
One is short-ranged, and one is of pion range |SH1 1211 UHl HOI IS2] • Each of them has an 
undetermined coefficient. In principle those coefficients should be fit to data (e.g. the 
deuteron magnetic moment, which is not exactly reproduced by the current J and the 
wave functions employed here) and the low-Q^ shape of the form factor. 

3 Unitary equivalence and the consistent treatment 
of 1/M corrections 

In this section we discuss the constraints imposed by Poincare invariance — or the low- 
energy manifestations thereof — on the Breit-frame isoscalar A^A^ charge operator. Recall 
that in HBxPT (pVo^^(q)lp) 

arises from a piece of C^^lf that has a fixed coefficient 
obtained via a Foldy-Wouthuysen transformation. Therefore this piece of the charge 
operator is a low-energy consequence of Lorentz covariance of (M'| J^|M). As such the 
contribution (^01) should be computed in a manner consistent with that used to derive 
the corrections to the one-pion exchange part of the A^A^ potential. Those 

corrections can be obtained from the chiral Lagrangian — specifically from the 1/M pieces 
m and the l/M^ pieces in C^^lf. But the relevant operators involve the energy of 
the individual nucleons, and so it is not immediately obvious how to convert them to 
contributions to an energy-independent quantum- mechanical potential. In fact, in the 
1970s and 1980s many techniques were developed by which quantum-mechanical operators 
could be obtained from a relativistic quantum field theory [HSl HH 112] • In all of these 
techniques there was freedom in choosing whether (and if so, which) nucleon lines to 
put on shell, as well as freedom in how to include meson retardation. As we shall see. 
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the choices made with respect to these two issues have an impact on the form of the 
operators (both V and Jo) that are obtained. Ultimately though, as long as operators 
and potentials are derived in a consistent way, the different choices are related by unitary 
transformations that leave matrix elements unaffected [HI |321 BH] • 

That unitary transformation is labeled by two parameters: u, which parameterizes the 
energy, k^, of the exchanged pion via 

*5 = (l-2.)^^^:^, (24) 

where p' (p) is the length of the relative three-momentum vector of the A^A^ system after 
(before) the meson exchange; and f3, which denotes a choice for the change in the nucleons' 
energy after absorption of the pion ^] : 

(p[,-Po)i,2 = (l-2/?)^^. (25) 

Note that in quantum mechanics energy is not conserved at each vertex, and so (Pq— Po)i = 
^0 = —{p'o ~ Po)2 need not necessarily hold. Indeed, it turns out that since we are in the 
A^A^ cm. frame the difference p'^ —po is the same for both nucleons once an energy shell is 
chosen. The full expression for the corrections to the one-pion-exchange potential in 

the case of arbitrary j3 and u can be found in Ref . . The main result for our purposes 
here is that if /? = | then the potential takes the form: 

/r.'lT/(l-)lr.\ ,a,a ^1 ' (P^ " P)^2 ' (p^ - p) f, p'^ + p" , ^ ( ^ \ \ 

^P'^ 'P^ = ~"^"^4g (p'-p)^ + m^ i'-^M^ + ^Mj- (26) 

This is the one-pion-exchange potential used in the N'^LO computation of Ref. ^Hl- (Cor- 
rections to "leading" two-pion exchange diagrams that are suppressed by 1/M are also 
included there, but are associated with pieces of the Jo which are of higher order than 
we work to here.) The computation of Ref. j2ZI employed the form for (p'jV^'-^'^-' |p) that 
corresponds to /3 = 0. All other potentials we discuss here (including the NNLO and NLO 
ones used in Ref. jTH|) employed the non-relativistic form of OPE, i.e. the result (|^. but 
without the additional factor in the round brackets. Meanwhile, all of the potentials we 
have used neglect retardation, which means they have set z/ = | /cq = 0. 

Consistent reduction of the contributions to the deuteron charge operator then leads 
to a more general result for diagram Fig. ^b) than that given in Eq. ()2()|1 [.41] : 



(pVf (q)|p) 



'1 '2 



8/2M 



\. _ . ffi -qtis ■ (p-p^ + q/2) 
, m2 + (p - p' + q/2)2 



l-^ ^^.(p_p^ + q/2)a,.(p-p^ + q/2) g . (p - + q/2) 
2 K + (p-p' + q/2)2]2 +1 J 



(27) 
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In Eq. fl20j) we obtained the result for the 0{eP^) piece of the charge operator that 
corresponds to (3 = and z/ = 1, because the field-theoretic manipulations used to arrive 
at Eq. (j^UI) assume that the fields represent physical particles, i. e. they are on-shell. The 
result ()27|1 may be obtained from Eq. ()2()|1 by applying a unitary transformation [HI |32] : 

4'\(3,u) = U\(3,u)j!t\0,l)U{f3,u), (28) 

where the form of U can be found in the original papers. The same unitary transformation 
generates consistent corrections to the one-pion-exchange part of the NN potential: 

VopeW, = U\P, z/)\/ope(0, l)U{P, u), (29) 

including the form ()2(i|l if the choice P = ^, u = ^ is adopted. This is not consistent with 
the choice made in obtaining Eq. ()20j) in Ref. ^1] because the NN potential of Ref. ^Hl 
was computed using the Okubo formalism developed in Ref. |13j. 

In Ref. |18j this issue of the choice made for /3 and u does not arise until the N^LO 
potential is derived, because in that paper, and in the earlier Refs. Epelbaum 
et al. chose to count M ~ A^. In doing this they were adhering to Weinberg's original 
argument as to why it is the two-nucleon-irreducible kernel — and not the A^A^ amplitude 
itself — which admits a chiral expansion. A^A^ intermediate states introduce factors of M 
in the amplitude for loop graphs, and if M ~ then the rath iterate of the one-pion- 
exchange potential is the same order as one-pion exchange itself IT^ . However, in 
Ref. [21] the need to iterate the one-pion-exchange potential to all orders was established 
without any reference to counting M ~ A^, being justified instead by the singular, and 
attractive, nature of the A^A^ force (see also Ref. ^3)- Therefore, while it is true that 
corrections to the A^A^ potential which are suppressed by powers of 1/M are often smaller 
than, e.g. those arising from excitation of the Delta(1232), in discussing electron-deuteron 
scattering we will consider a regime in which q/M can be sizeable. Therefore we follow the 
original HBxPT counting and take M ~ A. As we shall see, this counting is supported 
by the fact that the contribution (j2(J|) plays a significant, but not excessive, role in the 
deuteron charge and quadrupole form factors. 

If we count M ~ A the dilemma presented by the inconsistency between V and Jq 
arises already at O(eP^). A way out of this dilemma is provided by Eqs. and (|^. 
They guarantee that we will obtain the same result (up to 0(p'^/M^) corrections) for 
deuteron matrix elements (M'|Jo|M), regardless of what choices for (3 and v we make 
when constructing the operators V and Jq from the chiral effective field theory, provided 
that we consistently include the 0{p'^ /M"^) pieces of the potential V and the O(eP^) 
pieces of the operator Jq- Therefore in order to be consistent with the calculation of the 
corrections to Vope in Ref. [TH] we must adopt the choice (3 = \ m. the formula 
pTjl for \ If we do this, and also make sure to calculate one-pion exchange according 
to the formula (pUj) . then our results for matrix elements of the deuteron charge operator 
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will incorporate the low-energy consequences of Lorentz invariance, up to corrections of 
0{p^/M^) (higher order than we consider here). Note that the CD-Bonn potential is a 
different case, since the use of a pseudoscalar ttNN coupling means that there we have 
(3 = 0. Therefore in that case, and only in that case, we have used the expression ()20|) 
for the first part of Jq , with no modification by the factor of | that must be present if 
V^^'^^ is constructed with 13 = \. 

This still leaves us with the issue of how the /M'^ corrections in the one-pion exchange 
potential (j^UI) and the •p? jM"^ corrections to the nucleon kinetic energy operator are to be 
accounted for in the calculations using the NNLO and NLO wave functions of Ref. ^H] (or 
included in calculations with the Nijm93 wave function of Ref. j2^ or the wave functions 
of Ref. [2H1)- The original calculations of these wave functions did not include such effects, 
but since we count M ~ A here, we need to include them in order to have a consistent 
calculation of Gc and Gq to O(eP^). 

Starting from the Kamada-Glockle transformation we show in Appendix 1X1 the 
major part of these effects can be included by making changes to the short-distance pieces 
of the potential, and using a slightly modified wave function ^/^ in the computation 
of Gci Gq, and Gm- That wave function is related to the original non-relativistic wave 
function ip obtained in Ref. [TH1 by 

- ilWW [ m + 71r' + P j '"^ (yW^ - ■ (30) 

The solution of the non-relativistic Schrodinger equation for the wave function ?/', followed 
by the use of the formula (jHUj) to obtain the solution of the relativistic Schrodinger equa- 
tion, is the method by which we incorporate effects for the NLO and NNLO wave 
functions of Ref. ^H] , the wave function of Ref. [2^] i and the wave functions of Ref. ^28] . 
The effects of using ip, rather than the wave function, ip, to calculate electron-deuteron 
observables increase with photon momentum transfer |q|, but are small over the entire 
range for which the xET predictions can be trusted. At |q| = 700 MeV for the Nijm93 
wave function they change Gc by 6.3%, Gq by 1.2%, and Gm by 2.0%. (The effect on Gc 
is proportionately larger because 700 MeV is quite close to the form factor minimum.) 

4 Short-distance A^A^ charge operators and Qd 

So far we have obtained the deuteron two-body charge operator up to 0{eP^), or next-to- 
next-to- leading order. This is the order up to which the calculation we present here is fully 
systematic. In this section we discuss the role of contributions that are nominally higher 
order, and identify one particular mechanism that apparently could generate significant 
effects at = 0. We are particularly interested in this operator because "the Qd problem" 
that is present in all modern potential models (see, e.g Ref. 1211113) persists in the xET. 
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The problem is that all such calculations under-predict the value (H)) for by about 2-3% 
when they use a charge operator that includes all effects up to NNLO in the xET. The 
remaining discrepancy is large compared to the expected size of higher-order effects. 
It is also large compared to other discrepancies between theory and experiment in the 
^Si-'^Di channel of the NN system. 

And the situation is actually worse than this, because at 0{eP^) — one order higher than 
we are considering here — there are two-meson-exchange contributions to the deuteron 
charge operator. One might hope that these provide the missing strength in the E2 
response of deuterium at = 0. These diagrams are presently being calculated for finite 
Q^, and will be incorporated in a future computation of the charge and quadrupole form 
factors jini- However, it is already known that they do not give a sizeable contribution 
to the deuteron quadrupole moment. Park et al. 140J computed their effect on Qd using 
the AV18 wave function [47J, and found: 

^Qf = -0.002 fml (31) 

Therefore these effects will not repair the discrepancy between the calculated Gq{0) and 
the experimental Qd- 

At the next order, O(eP^), there are additional two-pion-exchange contributions to the 
deuteron charge. However, short- distance two-body currents that contribute to (rf) and 
Qd are also present, and are depicted in Fig. life). Even though it is suppressed by five 
powers of P relative to the LO result, the latter operator can have a noticeable impact 
on the quadrupole moment of deuterium, since the numerical value of Qd corresponds to 
a distance that is small on the typical scale of deuteron physics ~ 2 fm. The operator 
is 123: 

(pVf (q)|p) = |e|(l + n ■ ^2)^ (^1 ■ q^2 • q - ^-^a, ■ a^j , (32) 

and is designed to be diagonal in two-body spin and isospin and contribute only for 5* = 1 
and T = states. In the case of deuterium it represents an E2 photon inducing a ^Si ^Si 
transition. Such a transition is possible because the photon interacts with the total spin 
of the two-nucleon system through the two-body operator ()32|). The two-nucleon operator 
()32|) will be induced when high-momentum modes in the NN system are integrated out to 
obtain the low-momentum effective theory. It could also be induced when heavy mesons 
which can couple to a quadrupole photon in the requisite way are integrated out of the 
xET. This heavy-meson origin for the operator leads us to anticipate that with a scale Aq 
of about 1.2 GeV the coupling /14 will be of order 1. In particular, if we used resonance 
saturation in the A^A^ system ^H] to estimate the size of this operator the first mesonic 
current that would contribute to the operator would be the pai7 current pS] . 

We now give arguments which demonstrate that physics at roughly this scale could 
indeed remedy the discrepancy between the experimental Qd and the result found from 
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the mechanisms already discussed. Let us take the accepted number from "high-quahty" 
potential models = 0.270 fm^ (see, e.g. Ref. |lZj). Calculations with xPT wave 
functions obtain similar, or even slightly smaller, numbers pTl I2H1 1^- Then, we 
adopt AQ^a^ = 0.008 fm^ as an estimate for the NLO and NNLO corrections (which come 
mainly from the two-body operator Jq^"*). This leaves us with a remaining discrepancy 
between theory and experiment of 0.008 fm^, or about 3%. Inserting the operator (j32|) 
between deuteron wave functions we obtain its contribution to Qd as: 

^QT = ^-^\m\^ (33) 

where ip{0) is the deuteron wave function at the origin. While ipi^) is not an observable, 
neither is the dimensionless number h^. it is a wave-function- regularization dependent 
coefficient. Since only S-waves contribute at r = if we assume Aq = 1.2 GeV we can 
constrain the combination of ip{0) and to be: 



hi 



hm^ 

r—^O r 



6.5 fm~^ (34) 



where u{r) is the ^Si radial wave function. Therefore the operator ()32j) can be associated 
with physics at scales of about 1 GeV and still remedy the discrepancy between the 
theoretical value of Qd (including the meson-exchange contribution to the charge ipTjl ) 
and the experimental value (P). 

This higher-order effect has an importance greater than one would anticipate in Wein- 
berg's counting (|TT)|l not because it is unnaturally large, but because, from the point of 
view of that counting, Qd is unnaturally small [28 . This is hardly a surprise, since, at 
both leading and next-to-leading order, Qd is generated by one-body operators that con- 
nect the deuteron S-state to the deuteron D-state. Such effects are suppressed by the 
ratio Tj = An/ As = 0.0253(2) [50^ In contrast, the operator (jH^ is not ?7-suppressed 
and so its contribution to Qd is significantly larger than the naive estimate of ~ 0.1% 
leads us to anticipate. In this context it is worth noting that such an estimate for the 
short-distance effects in (r^) is completely validated by calculation The leading-order 
piece of (r^) is of the expected size ~ ^/ml, and two-body contributions, beginning with 
effects from Jg^'* and continuing through the two-pion-exchange mechanisms of Jq^"* and 
the CO-photon short-distance operator of Jq , give contributions of the expected size, 
approximately 0.3%. 

While this is reassuring, the (relatively) large impact of Jq'^ on Gq{0) means that 
we must ask how well we know this operator. Its value at = can be fixed by the 
requirement that it repair the discrepancy between theory and experiment for Qd- But 
at this stage we know nothing about its Q^ dependence. For the purposes of this paper 
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we will assume that this operator arises from heavy-meson exchange, and so model its 
dependence by: 




The uncertainty in the operator is now encoded as uncertainty in the scale A of its 
momentum variation. We anticipate A > 1.2 GeV, because there are no meson resonances 
below 1.2 GeV which, when integrated out of the theory, will yield this operator. The 
only danger with this reasoning is that two-pion-range mechanisms that occur at O(eP^) 
may ultimately prove to be responsible for the Qd discrepancy. This possibility is under 
investigation jlHl- However, evaluation of relevant processes in models which calculate, e.g. 
the role of AA components in deuterium, suggest that the dominant two-pion-exchange 
contributions to Jq'^ are not large enough to remedy the Qd discrepancy [HH E2] • 

As for an upper bound on the value A; the effects of the operator persist to higher 
as the scale of its momentum variation is increased. At the Q^'s considered here, the 
impact of this operator on observables when we choose A = 2 GeV is within 30% of what 
one would obtain at A = oo, so we will vary A between 1.2 and 2 GeV in order to assess 
the theoretical uncertainty of our calculation. We will see below that even with this range 
of variation our ignorance as to the precise value of A (or the precise function of that 
modulates the current Jq'^) is the dominant contribution to our theoretical uncertainty 
in the ratio Gc/Gq. 

Our goal in introducing the Q^-dependence (IS^ into our calculation is to assess the 
potential impact on our xET calculation from physics that is not explicitly included in 
it. The Q^-dependence of Jg'^^ will be modified by these sorts of effects, as well as by 
higher-order loop contributions that can be calculated in the xET. However, once such 
higher-order calculations are carried out the Q^-dependence of Jg^'' can presumably be 
constrained by input from electro-production in the single-nucleon sector. Therefore here 
we take the operator Jq as given. When we quote ranges for its impact on observables 
those ranges arise from the fact that {M'\J^^\m) varies when evaluated with different 
deuteron wave functions. 



5 Results for Gc and Gq 

In this section we present results for the matrix elements of the deuteron charge operator: 
Gc and Gq. 

In Figure|21we show the results for Gc when the NNLO [0(eP^)] operator for Jq is used 
(with nucleon structure included via Eq. (j22I))- The constants employed in evaluating the 
operator were qa = 1.29, U = 93.0 MeV, = 139 MeV, and M = 938.9 MeV. The 
dashed, dot-dashed, and solid lines show the range of predictions generated using NLO, 
NNLO, and N^LO wave functions with different regulator masses A and A. A list of the 
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Wave function 


Order 


A (MeV) 


A (MeV) 


001 


NLO 


400 


500 


002 


NLO 


550 


500 


003 


NLO 


550 


600 


004 


NLO 


400 


700 


005 


NLO 


550 


700 


101 


NNLO 


450 


500 


102 


NNLO 


600 


500 


103 


NNLO 


550 


600 


104 


NNLO 


450 


700 


105 


NNLO 


600 


700 


221 


N^LO 


450 


500 


222 


N^LO 


600 


600 


223 


N^LO 


550 


600 


224 


N^LO 


450 


700 


225 


N^LO 


600 


700 



Table 1: Values of the SFR cutoff A and the Lippmann-Schwinger equation cutoff A 
that are employed in the different wave functions of Ref. [TH] that are used to compute 
deuteron form factors in this work. The wave functions are in groups of five: first those 
generated with the NLO xET potential, then NNLO, then N^LO. 
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values of A and A that are chosen for the wave functions used here is given in Table ^ 
(which is adapted from Ref. [H])- At a given order in the expansion for the chiral potential 
these wave functions all have the same long-distance part, but the different scales at which 
spectral-function regularization is applied to obtain the A^A^ potential, and at which the 
exp(— p^/A^) regulator is applied to the potential before its insertion into the Lippmann- 
Schwinger equation, mean that they differ in their short-distance physics. Therefore the 
amount by which predictions for Gc vary once the order of the wave-function calculation 
is fixed gives us a lower bound for the impact of short-distance physics on our calculation. 




q (MeV) 

Figure 2: Results for the charge form factor Gc as a function of |q| = a/Q^. The dashed 
lines show the largest and smallest form factors obtained with the NLO wave functions of 
Ref. p^. The range of predictions with these wave functions is given by the horizontally 
shaded region. Similarly for the dot-dashed lines and the diagonally shaded region at 
NNLO, and the solid lines and the vertically shaded region at N^LO. Other results shown 
are for a wave function from Ref. with R = 1.5 fm (short-dashed line), and the 
Nijm93 and CD-Bonn deuteron wave functions (solid red and blue lines respectively). The 
experimental data is taken from the extraction of Ref. [22] and from Ref. upward 
triangles represent data from the T20 measurement of Ref. [Slj, open circle solid circle 
jSni, open squares |HZ], downward triangles j2H], rightward triangles (3^1, star jnOI, solid 
squares [HI], solid diamonds 

It is worth noting that the NLO potential only includes O(P^) corrections to V, and so 
the calculations labeled "NLO" in Fig.|21are limited by the accuracy of the A^A^ potential. 
The calculations labeled "NNLO" have O(P^) corrections included in the potential: the 
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same level of accuracy to which the operator Jq has been computed. In this respect only 
the computation with the NNLO wave functions is one that is carried out to a consistent 
order in both the wave functions and operators obtained using xET. The results with 
the NLO and N'^LO wave functions are shown for comparison. For the same reason we 
show results with the Nijm93 wave function, the CD-Bonn wave function, and the one- 
pion-exchange plus square well wave functions of Ref . [2H] . In each case consistent choices 
for (3 and v (as explained in Sec. E]) were employed. In the case of all but the CD-Bonn 
and xET N^LO calculations the impact of the p^/M^ pieces of the one-pion-exchange 
potential and the kinetic-energy operator has been assessed via the techniques described 
in Appendix Therefore all of these matrix elements have the same one-pion-exchange 
physics. A difference in their predictions is then either due to physics at the range of two- 
pion-exchange, or to physics at distances less than the scale where the chiral expansion 
can be used to reliably calculate the NN potential. 

In fact, Fig.Elshows us that all of these wave functions predict very similar form factors 
for |q| < 600 MeV. The most noticeable difference occurs around the zero of Gq — a region 
where, by definition, higher-order contributions cancel with lower-order contributions, and 
the calculation is therefore sensitive to the addition of higher-order effects. Most wave 
functions also produce a Gc in good agreement with the data compilation of Abbott et 
al. inn] for |q| < 600 MeV. An exception is the predictions using the N^LO wave function 
of Ref. PHI, which diverge from those of the other wave functions considered here at 
significantly lower . Note also that the N'^LO predictions seem to be significantly more 
sensitive to short- distance physics than is the case for the wave functions computed with 
NLO or NNLO chiral potentials, or even than the wave functions of Ref. f^, where only 
the result with the square well and surface delta function with = L5 fm is shown, but 
changing R to 2.5 fm produces a barely discernible change in the short-dashed line. 

It is possible that the situation with the predictions from the N^LO potential will 
improve when the pieces of Jq of O(eP^) and O(eP^) which are not included in this 
calculation are added to the current operator. But, even if this is the case, sizeable can- 
cellations between lower and higher-order effects are necessary if the N^LO wave function 
is to be used to describe electron-deuteron data at momentum transfers > 0.2 GeV^. 

One might argue that one does not expect the xET to work beyond this scale anyway. 
However, the chiral expansion developed here and in Refs. j2Hl 1121 1^ for the deuteron 
current operator still converges well at ~ 0.3 GeV^. In part this is because the impulse 
(leading-order) result for Gc can be written as: 



where jo is a spherical Bessel function and w{r) is the '^Di radial wave function. This 
means that — at least for the impulse-approximation piece of the matrix element — the 




(36) 
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relevant momentum scale at which the deuteron wave function is probed is not, in fact, 
|q|, but |q|/2 — half of the momentum transfer is taken away by the center-of-mass degree 
of freedom. In this context the failure of the N^LO wave functions' form-factor predictions 
to agree with the data when |q|/2 ^ 200 MeV is rather disturbing. 

In Fig. ini we show the results for Gq. As mentioned above, the shape produced by 
all wave functions is remarkably similar — a point which was exploited in an extraction 
of G^E^ in Ref. [EE]. This can be understood from the presence of j2, instead of the jo 
of Eq. in the co-ordinate space integral that gives the leading-order contribution to 
Gq. The pattern of convergence for Gq predictions with the order of the xET potential is 
interesting. The NLO band is quite wide, and its centroid is below the data. The NNLO 
band is very narrow, and its centroid agrees well with data out to 800 MeV. We stress that 
this is the consistent order for computation of the potential, given the current operator 
we have used. The N^LO band is then as wide, and below, the NLO band. In this plot 
predictions using the Nijm93 and CD-Bonn potentials are not shown. But they lie within 
the diagonally-shaded band that represents the range of predictions of the NNLO xET 
potential. 




q (MeV) 

Figure 3: Results for the charge form factor Gq (in units of fm^) as a function of |q| = 
\fQ^. Legend as in Fig. El apart from the absence of curves for the Nijm93 and CD-Bonn 
wave functions. These curves fall within the range of the dot-dashed lines, i.e. the results 
with NNLO xET wave functions. 

In order to remove the rapid fall-off in the plots of Fig. El and |21 and also provide a 
result for the ratio of form factors which will be measured at BLAST at the momentum 
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Experiment 


NLO:001-003 


NNLO:104-102 


N3L0:221-225 


Nijm93 


Qd (fm2) 0.2859(3) 0.278-0.280 


0.279-0.282 


0.270-0.274 


0.276 



Table 2: Deuteron quadrupole moment computed with our NNLO charge operator and 
different wave functions. Results are accurate to the number of digits shown. The ranges 
are generated by considering various values of A and A at a given order in the expansion 
for the xET A^A^ potential. The "extremal" wave functions are indicated in the top line 
of the table. 

transfers indicated by the asterisks, we show predictions for the ratio Gq/Gq in Fig. |3J 
These predictions are compared to data from the compilation of Ref. jSSl; as well as the 
more recent data set from Novosibirsk ^59^ ^. As can easily be gleaned from Fig. |3J the 
different wave functions considered give predictions that disagree at about the 5% level 
at = 0, i.e. they give different numbers for the deuteron quadrupole moment Q^- 

Numerical results for Q^, computed with the NNLO operator that was used to generate 
the predictions of Fig. 01 are presented in Table |21 Note that the variation in the results 
with the xET NLO and NNLO wave functions as the cutoffs A and A are varied is of the 
same size as the discrepancy between those predictions and the experimental value IQ. 

The fact that short- distance physics can affect the value of Qd at the 2-3% level needed 
to restore agreement between theory and data encourages us to incorporate the operator 
(see Eq. in our calculation. In doing so we adopt the Q^-dependence of Eq. (I35|) 
with A = 1500 MeV. For each wave function the value of is adjusted to yield the 
experimental value of Qd- The results of this procedure are presented in Fig. The 
NNLO, N'^LO, Nijm93, CD-Bonn, and one-pion-exchange-tail potential curves are only 
shown out to the momentum transfer where the contribution of the operator (j32p makes 
up 20% of the contributions from the preceding orders. This gives a way to estimate 
where the calculations with various wave functions become unreliable: they are unreliable 
when Jq'^ is no longer a small piece of the total Gq/Gq. Under this criterion most of the 
wave functions can give reliable predictions to |q| = 500-600 MeV, and below this value 
the wave-function dependence is markedly reduced as compared to what is seen in Fig. |3J 
Note that the lines at NLO are shown only to give an idea of the uncertainty coming from 
sensitivity to the choice of (A, A). Their predictions with this wave function for Gc/Gq 
beyond 600 MeV are provided only for recreational purposes. 

Lastly we focus on the region where the xET is reliable: |q| < 600 MeV. An enhanced 
view of this region is shown in Fig. El For each choice of A^A^ potential two curves are 

^The experimental error bars in this plot, and in the plots of Gc /Gq below, were obtained by summing 
the relative errors for Gc and Gq given in Refs. [S31 EH] • Some of the measurements of t^o and T20 from 
which these data came were quite precise, and so this procedure may well overestimate the size of their 
errors. 
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Figure 4: Results for the ratio Gc/Gq. As in Fig. |21the dashed hues and the horizontally 
shaded region show the range of results obtained with the NLO wave functions of Ref. [TB] . 
Likewise for the dot-dashed lines and the diagonally shaded region at NNLO, and the 
solid lines and the vertically shaded region at N^LO. Other results shown are for wave 
functions from Ref. 28 with R = 1.5 fm and R = 2.5 fm (short-dashed lines), and the 
Nijm93 and CD-Bonn deuteron wave functions (solid red and blue lines respectively). 
Upward triangles are data from the T20 measurement of Ref. [53], open circle solid 
circle jSH], open squares j^ZI, downward triangles [HH], rightward triangles jSH], star 60 , 
solid squares jHl], solid diamonds [H2j. The asterisks indicate the points where BLAST 
will extract this ratio from their t20 data. 




Figure 5: Results for Gc/Gq, after Jq is inserted with a coefficient adjusted to reproduce 
the experimental value of Qd and A = 1.5 GeV. Legend as in Fig. |3 Each curve is shown 
only up to the point where the Jq contribution is so large that the calculation is no 
longer reliable (with the exception of NLO). 
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q (MeV) 

Figure 6: Results for Gc/Gq, showing the variation that results from ignorance of the Q^- 
dependence of the operator Jq'^ (A = 1.2-2 GeV). Calculations are shown for NLO (long 
dashed), NNLO (dot-dashed), N^LO (solid green), Nijm93 (solid red), CD-Bonn (sohd 
blue) and R = 1.5 fm square well + one-pion exchange (short dashed) wave functions. 
The vertical lines indicate a reasonable range for the theoretical prediction at each of the 
points where BLAST will have data. The experimental data is from Refs. [S3] (upward 
triangle), (open circle), [3^] (solid circle), and jST] (open square). Note change of 
scale as compared to Fig. 



|q| (MeV) 


Gc/Gq (fm2) 


Error bar (fm^) 


368.4 


3.11 


0.11 


403.9 


2.99 


0.13 


439.3 


2.86 


0.16 


474.8 


2.71 


0.18 


514.2 


2.53 


0.22 


559.5 


2.29 


0.26 


606.8 


2.02 


0.30 



Table 3: Predictions for the ratio Gc/Gq at the values of |q| where this quantity will 
be measured at BLAST. The error bar and central values displayed here are obtained via 
the procedure discussed in the text. 
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shown: the upper one of the pair corresponds to choosing A = 1.2 GeV when evaluating 
the operator Jg and the lower one corresponds to choosing A = 2 GeV. A conservative 
estimate for the impact of short-distance physics which is not well-constrained in this 
xET calculation is given by combining the uncertainties from (A, A) variation and the 
uncertainty coming from lack of knowledge about the momentum dependence of Jq • The 
black bars then represent the range of possible values that the xET predicts for Gq/Gq 
at the kinematics where there will be data from BLAST. These ranges are reproduced 
in Table El The error is about ±3% at the lowest BLAST point and increases with Q^, 
as it should. Note that we have not included the N^LO xET wave function, or the NLO 
xET wave function, in generating these predictions. As already discussed, the predictions 
for Gc and Gq with the N^LO wave function deviate already from the extant data at 
quite low Q^, while the NLO wave function is, in the xET sense, less accurate than 
the operator being used here. The predictions obtained with these wave functions are, 
however, within 2cr, if the theoretical error bars we have obtained are taken to have the 
usual one-standard-deviation interpretation. 



6 Results for Gc/Gm 

BLAST will also measure the ratio Gq/Gm- Predictions for that observable are provided 
in Fig. 13 We do not show any experimental data in Fig. [3 because, as far as we can glean 
from the literature, all previous data on Gm comes from different data sets to that used 
for the extraction of Gq and Gc [SB]- Therefore in general data on Gc and Gm are not at 
the same and have different systematic errors. The BLAST data set will be pioneering 
in this regard. 

The calculations displayed in Fig. [3 are accurate to relative order P^, although the 
Gc used here is actually computed up to relative order P^. Once again the shape of 
the low-momentum part of the ratio is fairly wave-function independent, but the value at 
= changes as we move through the different wave functions used for the computation 
of Fig. due to short-distance contributions to Hd being different for different wave 
functions. However, even without renormalization there is a robust prediction for the 
ratio out to ~ 0.1 GeV^. The robust prediction is that GcjGu is approximately flat. 
This would be exactly the case in the absence of relativistic, meson-exchange, and nucleon- 
structure contributions to the operator, and if wir^ = 0. The relativistic corrections at 
0{P^) are negligible at < 0.1 GeV^, and the meson-exchange piece of the charge 
operator is higher order than we are attempting to calculate the ratio Gc/Gm to. As far 
as the operator is concerned this leaves only the nucleon-structure effects, which depend 
on the ratio: 
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-- NLO:001-003 
■- NNLO:102-104 
— CD-Bonn 




q (MeV) 



Figure 7: Results for the ratio Gc/Gm- Calculations shown are for extremal NLO (long 
dashed), extremal NNLO (dot-dashed), extremal N^LO (solid green) (solid red), CD- 
Bonn (solid blue) and i? = 1.5 fm square well -|- one-pion exchange (short dashed) wave 
functions. 



If a strict chiral expansion is used for the form factors then this ratio depends (again, at 
this order) on the difference of the isoscalar magnetic and charge radii, and amounts to a 
< 10% effect at = 0.1 GeV^. Even though taking the ratio Gc/Gm does not (as it did 
in the case of the ratio of the previous section) eliminate nucleon-structure effects, it does 
reduce their size. Meanwhile, the effects of w grow with and so at < 0.1 GeV^ it 
is thus not particularly surprising that Gm/Gq is, to quite a good approximation, flat. 

Given the extent of the variation in the prediction for the ratio beyond |q| = 500 MeV 
it is difficult to believe that the NLO predictions for the ratio shown here are reliable 
beyond that point. This situation could improve once the NNLO pieces of the operator 
J were computed, but short-distance pieces of J appear already at that order. Therefore 
it is a prediction of the chiral expansion that this ratio will be more sensitive to short- 
distance physics than is Gc/Gq. The position of the minimum in Gm is known to be very 
sensitive to such short-distance physics IM | IM j I65j . In this context it is worth noting 
that the minimum for the N^LO wave functions is already at |q| ^ 800 MeV — much lower 
than in any of the calculations of Refs. IHSl lEll tEH] and indeed, much lower than the 
experimental data allows the position of the minimum to be 
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7 Conclusion 



We have used the xET isoscalar charge operator in the nucleon-nucleon sector computed 
to NNLO (including a consistent treatment of the 1/M pieces of the charge), and the wave 
functions of Ref. [TH], to obtain the form-factor ratios Gc/G^^\ Gq/G^^\ Gc/Gq. These 
ratios test xET's abihty to describe deuteron structure. We confirm and extend the finding 
of Ref. ISD, that the NLO and NNLO xET wave functions, combined with the NNLO 
Jo, yield results for these ratios that agree — within the experimental uncertainties — with 
the extractions of Ref. ^ for < 0.35 GeV^. In contrast, the N^LO wave function of 
Ref. fHj, when used in conjunction with the N^LO charge operator, produces form factors 
that depart from the data at Q'^ ^ 0.2 GeV^. 

In light of the upcoming release of data on ^20 from BLAST we examined the ratio 
Gc/Gq in detail. We found variation in the xET predictions for this ratio at = 0, 
and also found that — even allowing for this variation — xET is in disagreement with the 
experimental value for this quantity. This phenomenon — the "Qd problem" — is familiar in 
high-precision NN potential models with the modern value of the ttNN coupling constant. 
In xET its solution arises naturally through a higher-order two-body current that couples 
exclusively to quadrupole photons. We added this operator to our analysis, and showed 
that when we do so the xET predictions for Gc/Gq (with the NNLO wave functions) 
fall within a narrow band out to ~ 0.35 GeV^ (see Fig. We also performed the 
calculation with the same charge operator and potential-model wave functions that have 
a one-pion-exchange tail identical to that of the xET potential I2Z1 12H|- We found 
that these wave functions make the band obtained at NNLO in xET about a factor of 
two wider. We conservatively adopt the full width of that band as representative of the 
theoretical uncertainty in our calculation. 

Meanwhile, the xET predictions for Gc/Gm^ which will also be measured at BLAST, 
are not reliable to as high a, Q^. In saying this, it should, in fairness, be pointed out that 
J has not been computed to as high an order as Jq, and it could therefore be that Gc/Gm 
can also be well described to = 0.35 GeV^ in xET once O(eP^) corrections to J are 
included. This is a topic for future work. Another topic for future work is the inclusion 
of the O^eP"^) pieces of Jo that were already computed in Ref. jlHI at = in the 
finite- 1 q| calculation [1^. In addition, the operators and the xET wave functions used in 
this paper to make predictions for the BLAST data can be further tested by comparing 
their predictions with experimental results for deuteron electro-disintegration — although 
this will require the computation of the isovector pieces of the operators. 

Irrespective of such future efforts, one thing is already clear from Fig. |B1 When the 
theoretical predictions for Gc/Gq are renormalized in the manner we advocate here, the 
theoretical uncertainty in Gc/Gq for < 0.3 GeV^ is less than the uncertainty in the 
data. This makes the low-Q^ data from BLAST all the more crucial, since it will provide 
an important test of xET's ability to organize contributions to deuteron observables, and 
its ability to use that organization to provide estimates of the theoretical uncertainty in 
a given calculation. 
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A Including p^/M'^ corrections in "non-relativistic" 
wave functions 

One way to include relativistic corrections to the nucleon kinetic energy operator in the 
xET is to solve the "relativistic Schrodinger equation" ^H] 



Vp^Tm^ - mJ ri'{p) + j ^ v,y{p,p')rAp') = E^t\p), (38) 

with the (partial- wave decomposed) potential Vii?{p,p') the one that is obtained from 
the xET using the counting rules explained in Section [T] To facilitate computation of 
the p^/M^ (and beyond) corrections to the nucleon kinetic energy operator we adopt the 
Kamada-Glockle transformation and define a new relative momentum p such that 



p2 
2M 

The inverse transformation is then: 



M2 + p2 _ M. (39) 



P = P\I^ + ^^=P9{P)- (40) 

As shown by Kamada and Glockle, any "relativistic" Schrodinger equation that employs 
the kinetic energy operator with the relativistic form on the right-hand side of Eq. (j39|l 
can be recast as a non-relativistic Schrodinger equation, i.e. jlSl HE] 

^'-M'{P) + / T^C yuKP.PW{p') = E^!'ip), (41) 



2M^' J (2n) 
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where the potential V^f/ is obtained from VJ^/ via: 

Vinp,p') = Hp)Viy{pgiP),P'g{p'))h{p'), (42) 

and 

^(p) = h{p)-^{pg{p)). (43) 

Here the factor h is introduced to ensure that, as long as ip is normalized to one, ip is also 
normalized to one. Calculation of the Jacobian associated with the transformation ()39|) 
yields: 



A(p) 



^^'1 + ^)9®- (44) 

We can also work this procedure in reverse, i.e. interpret solutions of the non-relativistic 
Schrodinger equation as solutions to a dynamical equation with a relativistic kinetic energy 
operator and a relativistic potential. (Something similar is done to obtain the A^A^ cm. 
frame Hamiltonian in approaches to electron-deuteron scattering based on Relativistic 
Hamiltonian Dynamics, see, for example, Ref. [HZI-)- Here we will make this interpretation 
for the deuteron wave functions obtained with the NLO and NNLO chiral potentials. 

Suppose that V is the chiral potential computed (at NLO or NNLO) without p^ / M"^ 
corrections, i.e. 

V^iKP.P') = f{P)M',i,{p,p')f{p'), (45) 

where M. is the (partial- wave decomposed) sum of A^A^-irreducible diagrams that is 
computed at NLO and NNLO using spectral-function regularization in Ref. PH| and 
f{p) = exp(— p^/A^) is the "Lippmann-Schwinger equation regulator" used there. A 
potential V that includes p'^/M'^ corrections and is associated with the same sum of Feyn- 
man diagrams can be constructed by inverting Eq. (jl^ . We find that this potential, 
which when inserted in the relativistic Schrodinger equation (j38p will be phase-equivalent 
to V of Eq. (USD, is: 



ViAp,p') = jI^^^u' [\l2M^M^+p^-2M^, ^J2M^ + p'^ ^2Apj (46) 

We have not bothered to distinguish between p and p in the functions / and h, because 
we now define a new regulator function f = f /h to obtain: 



Vii^{p,p') = f{p)M"l, l^pM^M^ + p2 _ 2M2, pM^M^ + p'^ - 2M^ j f{p'), (47) 

up to terms that are of order p'^/M'^. Since we count M as the same order as the cutoff 
scale A absorbing the functions h into the definition of the regulator in this way makes 
physical sense. Apart from these short-distance effects the only difference between the 
potentials V and V is the use of the "stretched" momenta in V. 
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If we wish we can also incorporate the factor (l — that distinguishes Eq. 

from the non-relativistic treatment of one-pion exchange used at NLO and NNLO in 
Ref. [IB]. To do this we redefine the regulator again, writing 

V,P(P,P') = 7(p) fl - My, U2M^K-P+pi-2M\ ,J 2M ,J + p'^ - 2mA 

where, up to the order to which we work: 



1 - /(P'),(«) 



If one ignores the regulator /, then the factors in front of the one-pion-exchange part of 
the potential in Eq. ()48|) agree with Eq. (pUj) and the expression for the two-pion-exchange 
pieces has been modified only at N^LO and beyond. As for the short-distance pieces of 
the potential, the changes in the way the chiral NN potential is regularized [f ^ f ^ f) 
mean that the A^A^ LEG C2 will have to be re-defined: its value is shifted by a term of 
0(1/M2). However, C2 is fitted to data, so no change in observables should result from 
doing this. 

Therefore we have shown that, via modifications of the (unobservable) short- distance 
part of the potential that do not affect the low-energy observables, we can absorb much 
of the "relativistic effects" that would enter the NLO and NNLO chiral potentials were 
we to count M ~ A. The only such effect that cannot be absorbed by a redefinition of 
the regulator used when solving the Lippmann-Schwinger equation with the potential V 
is the different arguments at which M. is evaluated in Eq. pH|l as compared to Eq. P^jl . 
This difference is an NLO effect for the part of V coming from one-pion exchange, where it 
represents the interplay of the relativistic kinetic energy operator and the (leading-order) 
one-pion-exchange potential. This interplay is not included in our calculations, but all 
other relativistic effects are, and as reported in Sec. El they prove to be small. 

The main difference between interpreting the expression ()45|1 as a non-relativistic po- 
tential and interpreting it as resulting from a relativistic calculation after application of 
the Kamada-Glockle transformation therefore resides in the momentum arguments that 
enter the Schrodinger equation. In the first case the wave function that is the solution 
of that equation will be a function of p. In the second case it will be a function of 



p = Y 2M a/M^~+~p^ — 2M2. The relationship between the two wave functions is given by 
Eq. (jini). This can be inverted to yield Eq. (jHUj) . which we rewrite here: 




^*' = V7JP7F mTtS^I ^.|V2MVM^ + P-2MM. (50) 
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Thus, if we are provided with the wave-function ip that is obtained by solving the non- 
relativistic Schrodinger equation with the potential V, then we may derive from that the 
wave function ip, that is the solution of the relativistic Schrodinger equation p8|l with 
a phase-equivalent relativistic potential V fj48|) . with V constructed to include all the 
corrections to one-pion-exchange that must be present in the NLO and NNLO potentials 
if we count M ~ A. This procedure is not necessary in the case of the N^LO potential. In 
that case the corrections which are connected with the current were explicitly 

computed in Ref. [THj . 
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